Charmonium in a hot, dense medium 



David Blaschke^'"^ 

^Institute for Theoretical Physics, University of Wroclaw, 50-204 Wroclaw, Poland 
^Bogoliubov Laboratory for Theoretical Physics, JINR, 141980 Dubna, Russia 

DOI: http://dx.doi.org/10.3204/DESY-PROC-2009-07/blaschkejdavid 

In this lecture we apply a thermodynamic Green function formalism developed in the 
context of nonrelativistic plasma physics for the case of heavy quarkonia states in strongly 
correlated quark matter. Besides the traditional explanation of charmonium suppresion 
by Debye screening of the strong interaction, we discuss further effects of relevance when 
heavy quarkonia states propagate in a medium where strong correlations persist in the 
form of hadronic resonances. These effects may be absorbed in the definition of a plasma 
Hamiltonian, which was the main result of this work. This plasma Hamiltonian governs 
the in-medium modification of the bound state energy levels as well as the lowering of the 
continuum edge which leads not only to the traditional Mott effect for the dissociation 
of bound states in a plasma, but can also be applied for a consistent calculation of the 
in-medium modification of quarkonium dissociation rates. 

1 Introduction 

In developing a theoretical approach to heavy quarkonia as messengers of the deconfinement / 
hadronization transition of a quark-gluon plasma formed in a heavy-ion collision, we should 
aim at a unifying description where hadrons appear as bound states (clusters) of quarks and 
gluons. The situation is analoguous to the problem of two-particle states in QED plasmas 
where a well-developed theory in the framework of the Green function technique exists. These 
methods have been widely developed for the case of the hydrogen plasma, where the electrons 
and protons as the elementary constituents can form hydrogen atoms as bound states of the 
attractive Coulomb interaction. The problem is tractable analytically for the isolated two- 
particle system, with a discrete energy spectrum of bound states and a continuous spectrum of 
scattering states. Higher complexes, such as molecular hydrogen can also be formed. 

In a many-particle system, the problem of bound state formation needs to account for 
medium effects. They give contributions to a plasma Hamiltonian 

^pl _ ^Hartrco _^ ^Fock ^ ^Pauli _^ jjMW ^ ^Dcbyo _^ jjpp ^ jjvdW _j_ 

where the first three effects, the Hartrec- and Fock energies of one-particle states and the 
Pauli blocking for the two-particle states, are of first order with respect to the interaction and 
determine the mean-field approximation. The following two terms of the plasma Hamiltonian 
are the MontroU-Ward term giving the dynamical screening of the interaction in the self-energy, 
and the dynamical screening (Debye) of the interaction between the bound particles. These 
contributions are related to the polarization function and are of particular interest for plasmas 
due to the long-range character of the Coulomb interaction. In a consistent description, both 
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terms should be treated simultaneously. The last two contributions to the plasma Hamiltonian 
are of second order with respect to the fugacity: the polarization potential, describing the 
interaction of a bound state with free charge carriers, and the van der Waals interaction, 
accounting for the influence of correlations (including bound states) in the medium on the 
two-particle system under consideration, see [1] [2] . 

Approximations to medium effects in the self-energy and the effective interaction kernel have 
to be made in a consistent way. resulting in predictions for the modification of one-particle and 
two-particle states. On this basis, the kinetics of bound state formation and breakup pro- 
cesses can be described, establishing the ionization equilibrium under given thermodynamical 
conditions [3]. Coulomb systems similar to the hydrogen plasma arc electron- hole plasmas in 
semiconductors [4], where excitons and biexitons play the role of the atoms and molecules. 
Other systems which have been widely studied are expanded fluids like alkali plasmas or no- 
ble gas plasmas, see [1] and references therein. Applications of the plasma physics concepts 
for cluster formation and Mott effect to the rather short-ranged strong interactions have been 
given, e.g., in [51 [5] for nuclear matter and in [71 [5] for quark matter. 

In this subsection, we want to discuss basic insights from these investigations of bound state 
formation in plasmas, as far as they can concern our discussion of heavy quarkonia formation 
in hot and dense matter. Before going more into the details, let us mention them. Bound state 
properties remain rather inert to changes of the medium since the self-energy and interaction 
effects partially compensate each other in lowest order of density. Also, the smaller size of the 
bound states matters in this respect. The compensation does not hold for continuum states, 
being influenced by self-energy effects only, so that a lowering of the in-medium ionization 
threshold must occur which leads to a strong enhancement of the rate coefficients for bound- 
free transitions and to a sequential "melting" of different bound state excitation levels into the 
continuum of scattering states at corresponding critical plasma parameters (Mott effect [5]), 
until even the ground state becomes unbound. 

The theory of strongly coupled plasmas has been developed also for strong nonideality, 
where the formation of clusters in the medium need to be taken into account. This situation 
is similar to that of a hadronizing quark-giuon plasma and we will therefore refer to cluster 
expansion techniques as the theoretical basis. 

2 Bethe-Salpeter equation and plasma Hamiltonian 

The most systematic approach to the description of bound states in plasmas uses the Bethe- 
Salpeter equation (BSE) for the thermodynamic (Matsubara-) two-particle Green function 

Gab = Cj^t + Kab Gab = + G"^ Tab G°f, , (2) 

which is equivalent to the use of the two-particle T-matrix Tab and has to be solved in conjunc- 
tion with the Dyson equation for the full one-particle Green function, 

Ga = GJ] + G'^a'^aGa , (3) 

defined by the dynamical self-energy Sa(?', w) and the free one-particle Green function Gq(p, uj) = 
[lo — ea(p)]~^ for a particle of species a with the dispersion relation ea{p) = y^P^ + iTia ~ 
ma+py{2ma), see Fig. [H 

The BSE contains all information about the spectrum of two-particle bound as well as 
scattering states in the plasma. A proper formulation of the plasma effects on the two-particle 
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Figure 1: The two-particle problem in the medium. Dyson equation (left) and Bethe-Salpeter 
equation (right) need to be solved in consistent (conserving) approximations for self-energy (E) 
and interaction kernel {K). 

spectrum is essential to understand why bound and scattering states are influenced in a different 
way by the surrounding medium, leading to the Mott-effect for bound states. We give here the 
essence of a detailed discussion to be found in Ref . [2] . 

The homogeneous BSE associated with ^ can be given the form of an effective Schrodinger 
equation for the wave function V'a6(pi,P2, z) of two-particle states in the medium [3] 

^{[£a{Pl) +£b{P2) - z]Sq.o - Vab{q)} IpabiPl +q,P2 ~ q, z) = 
Q 

= ^H^l{.Pi,P2,q,z)ipab{pi + q,P2 - q,z), (4) 

Q 

where a, b denote a pair of particles with 3-momenta pi and P2 which transfer a 3-momentum 
2q in their free-space interaction Vabiq) and ^ is a complex two-particle energy variable. The in- 
mcdium effects described by ([2|) have been singled out in the definition of a plasma Hamiltonian, 
containing all modifications beyond the two-body problem in free space [21 14] 

H^iipuP2.q,z) ^ Va,{q) [Nab{puP2) - 1] ^y^Vatiq) [Nab{pi^q\p2-q) - l]S,^o 

^ V ' , 

(i) Pauli blocking ^ ^ ^ 

(ii) Exchange self— energy 

-I- ^Vab{j>\-,P2-,q-,z)Nab{j>\-,P2) - A Kb (Pl , P2 , , ^)^a6 (Pl + , P2 " ?') f^g.O ■ (5) 
' ^ ' 9^ 

(iii) Dynamically screened potential v ^ y 

(iv) Dynamical self— energy 

Here, AT4h(pi,p2j 9, z) = Kab{pi,P2, q, z) — Vabisi) stands for the in-medium modification of the 
bare interaction potential to a dynamically screened interaction kernel Kab{pi,P2,q, z). The 
effects of phase space occupation are encoded in the function Nab{pi,P2), which for the case 
of an uncorrelated fermionic medium takes the form of the Pauli blocking factor Nabipi,P2) = 
1 - faiPi) - fb{p2), where fa{p) = {exp[(ea(p) - f^a)/T] + 1}"^ is the Fermi distribution and 
jia the chemical potentential of the species a. Eq. ([?]) is a generalization of the two-particle 
Schrodinger equation, where on the left-hand side the isolated two-particle problem is described 
while many-body effects due to the surrounding medium are given on the right-hand side. The 
in-medium effects named in the plasma Hamiltonian ([T]) can be obtained from the one derived 
in the Bethe-Salpcter approach ([5|) upon proper choice of the interaction kernel Kab so that 
Eq. ll]) appears a sa special case of Eq. ([5|) . 

The influence of the plasma Hamiltonian on the spectrum of bound and scattering states 
can be qualitatively discussed in perturbation theory. Since bound states are localized in 
coordinate space, their momentum space wave functions extend over a finite range A and we 
may assume them to be g-independent: ipabipi + q,P2 ~ q, z = Enp) ~ ^ah(pi,P2, z = Enp) for 
small momentum transfer q < A and to vanish otherwise. Assuming further a flat momentum 
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dependence of the Pauli blocking factors Nab{pi + q,P2 — 9) ~ Nab{pi,P2) for small q where the 
interaction is strong, we obtain a cancellation of the Pauli blocking term (i) by the exchange 
self energy (ii) and of the dynamically screened potential (iii) by the dynamical self-energy 
(iv). Therefore, the bound states remain largely unmodified by medium effects. For scattering 
states which are extended in coordinate space and can be represented by a delta function in 
momentum space, the above cancellations do not apply and a shift of the two-particle continuum 
threshold results. For this mechanism to work it is important that approximation schemes for 
the sclf-encrgy and the interaction kernel have to be consistent as, e.g., in the conserving scheme 
of (f>-derivable theories [TUl . 

Summarizing the discussion of the plasma Hamiltonian: bound state energies remain un- 
shifted to lowest order in the charge carrier density while the threshold for the continuum of 
scattering states is lowered. The intersection points of bound state energies and continuum 
threshold define the Mott densities (and temperatures) for bound state dissociation. 

When applying this approach to heavy quarkonia in a medium where heavy quarks (either 
free or bound in heavy hadrons) are rare, then Nab = 1 so that both, (i) and (ii) can be safely 
neglected. The effects (iii) and (iv) stem from the dynamical coupling of the two-particle state 
to collective excitations (plasmons) in the medium. In the screened potential approximation, 
the interaction kernel is represented by V^ij{piP2, q, z) = V^^^{q, z)5 p_p^j^.p^52q,pi-p2 with 

yi((Z,z) = Vabiq) + Vab{q)\iab{q,z)Vi{q,z) = Vab{q)[l - ^abiq, z)Vab{q)r^ , (6) 

with the total momentum P and the momentum transfer 2q in the two-particle system. The 
most frequently used approximation for the here introduced polarization function Ilabiq,z), 
or for the equivalent dielectric function eabiq,z) = 1 — Habiq, z)Vabiq), the random phase 
approximation (RPA). In the next two paragraphs we discuss the static, long wavelength limit 
of the RPA and its generalization for a clustered medium. 

Example 1: statically screened Coulomb potential. The systematic account of the 
modification of the interaction potential between charged particles a and b by polarization of 
the medium is taken into account in the dynamical polarization function Ilab{q,z), which in 
RPA reads [2] 

d^P fa jE;) - fa{E;_^) 

(2^)3 - E-_^ -z ■ ^ ' 

For the Coulomb interaction, corresponding to the exchange of a massless vector boson, the po- 
tential is obtained from the longitudinal propagator in the Coulomb gauge is Vabiq) — ^a^b/q^- 
For a recent discussion in the context of heavy quark correlators and potentials see, e.g., pnil2| . 
Due to the large masses of the constitutents in the heavy quarkonium case, one may use a Born- 
Oppenheimer expansion to replace the dynamically screened interaction by its static (z = 0) 
and long-wavelength (q 0) limit. For nondegenerate systems the distribution functions are 
Boltzmann distributions and their difference can be expanded as 

fa{E;) - fa{E;_^) = e-^."/^ (1 - e-(^^.-^.")/^) « ~faiE;){E; - E^_^ - z)/T , (8) 

so that the energy denominator gets compensated and the polarization function becomes 

n^r(9,-) - -^'-f J ^Ue;) = -Sab^ . (9) 



Il^^^{q,z) = 2Sab I 
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The corresponding dielectric function ejj^ (5, oj) takes the form 



Hm£^P^(g,0) = l + 4> Ml, = i^eX(r). (10) 

a 

The screened Coulomb potential in this approximation is therefore V^f (q) = Vab{q) / 0) = 
^a&hl^'f' + /^l))- In this "classical" example of the statically screened Coulomb interaction, the 
contribution to the plasma Hamiltonian is real and in coordinate representation it is given by 

AK,(r) = -^(e-^°'- - 1) « a^Ln - %^lr , (11) 
r 2 

where a = e^/(47r) is the fine structure constant. For the change in the Hartree self-energy of 
one-particle states due to Debye screening we can perform an estimate in momentum space 
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This entails that to lowest order in the density the shift of the one-particle energies (continuum 
edge of unbound states) + St, = —a^jj compensates the contribution due to the screening 
of the interaction pT|) 

Aab « a^iD = 0{^na^ Q) , (13) 
in the wave equation (|4]). For the shift of the bound state energy levels follows [2l [T3] 

AE,,^^-^^ll{r),,, = 0{nal,) , (14) 

where aB,o = 1/(q! rn) is the Bohr radius. 

The Debye mass fj,D, equivalent to the inverse of the Debye radius td characterizing the 
effective range of the interaction, depends on the square root of the density n{T) of charge 
carriers. It is this different response of bound states and scattering continuum to an increase of 
density and temperature in the medium which leads to the Mott effect (see, e.g., Refs. in [3] and 
[T]) for electrons in an insulator: bound states of the Debye potential can only exist when the 
Debye radius is larger than Mott = 0.84 a^.o [II]- This entails that above a certain density 
even the ground state electrons become unbound and form a conduction band, resulting in an 
insulator-metal transition also called Mott-transition. Further details concerning this example 
can be found in Ref . [15] . 

In complete analogy to this electronic Mott effect it is expected that in hadronic matter 
under compression the hadrons as bound states of quarks undergo a Mott transition which 
results in a phase transition from the color insulating phase of hadronic matter to a color 
conducting or even color superconducting phase of deconfincd quark matter. This applies to 
light hadrons as well as to heavy quarkonia, whereby due to the different scales of Bohr radii 
the Mott dissociation of heavy quarkonia occurs at higher densities than for light hadrons. 

In most approaches the quark self energy effects arc neglected and one is left with the only 
medium effect due to a statically screened potential. This has the consequence that in such a 
picture the continuum edge of the scattering states remains unshifted and due to the lack of 
compensation the effective interaction leads to a strong medium dependence of the bound state 
energies (masses). For the electron- hole plasma in highly excited semiconductors it could be 
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shown experimentally, however, that the compensation picture is correct and the bound state 
energies remain almost unshifted |16j . 

One may of course absorb the self-energy effects into a redefinition of the effective interaction, 
by adding a homogeneous mean-field contribution. This is equivalent to the use of the Ecker- 
Weitzel potential |T7] 



It is interesting to note that recent investigations of the screening problem in the context 
of Debye-Hiickel theory [TH] and QQ correlators [TTl [T^] have obtained this continuum shift 
(— a^u) as a homogeneous background field contribution. According to the above lesson from 
plasma physics, however, this contribution should be attributed to the self-energy of the con- 
stituents rather than to the interaction kernel, since it determines the shift of the continuum 
edge. 

For the development of a comprehensive approach to heavy quarkonia in hadronizing hot, 
dense QCD matter another insight from plasma physics may be of relevance and will be dis- 
cussed next: the effect of strong correlations (bound states) in the medium. To this end, the 
bound states will be treated like a new species occuring in the system. Accordingly, addi- 
tional diagrams have to be taken into account which stem from a cluster expansion of the 
interaction kernel Kab and the corresponding self-energy Ea, see Figs. [2]|4l This leads in the 
plasma Hamiltonian H^^ to a generalization of the self-energy contributions (cluster-Hartree- 
Fock approximation), the distribution functions in the Pauli-blocking factors and the dynamical 
screening (cluster- RPA). The van-der-Waals interaction in Eq. ^ appears naturally as a con- 
tribution in the cluster expansion, describing polarization effects due to bound states in the 
medium. 

3 Cluster expansion for quarkonia in correlated medium 

In the vicinity of the plasma phase transition, correlations play an important role and their 
proper accounting requires rather sophisticated theoretical methods such as cluster expansion 
techniques. For the problem of charmonium in dense hadronic matter at the deconfinement 
transition, i.e. in the strong coupling case, wc suggest a systematic Born series expansion of 
collisions with free and bound states in the surrounding matter so that all terms linear in the 
density of free particles and bound states are taken into account. 

Wc describe the cluster expansion here in terms of its diagrammatic expressions for the 
interaction kernel and the corresponding self-energy. The 1"* Born approximation diagrams 
of this expansion are given in Fig. [2l see also the monograph [2] . The wavy line denotes the 
dynamically screened interaction V^^, which in a strongly correlated plasma receives contribu- 
tions from the polarization of the medium beyond the RPA, denoted as generalized (cluster-) 
RPA in Fig. HI see Ref. [Tn]. Bound and scattering states are described consistently in the 
two-particle T-matrices. For a generalization to higher n-particlc correlations, see [SJ I20j and 
the monograph [2] . The diagrams containing T-matrices do not contribute to the charmonium 
spectrum as long as the densities of the charmed quarks and of charmed hadrons in the medium 
are negligible. This is the situation expected for the rather low-energy CBM experiment. For 
the discussion of charmonium production at RHIC and at LHC the inclusion of these terms can 
be invoked. 

At the 2"'' Born order, we distinguish two classes of collisions with light clusters (hadrons) 
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Figure 2: Cluster expansion for interaction kernel for the two-particle problem in a strongly 
correlated medium (upper equation) and the corresponding self-energy (lower left equation) 
with a dipole ansatz for the vertex (lower right equation). 




Figure 3: Alternative way of drawing the diagrams for the cluster expansion of the interaction 
kernel and the corresponding self-energy of Fig. [2] in a form familiar in plasma physics and 
nuclear physics. 



that can give rise to spectral broadening of the charmonia. The first class concerns hadron 
impact without quark rearrangement inducing transitions to excited states, shown in the left 
panel of Fig. [5l These processes have been considered for charmonium-hadron interactions 
within the operator product expansion techniques following Pcskin and Bhanot [2T1 [22] , see [23^ 
124] . The result is a deformation of the charmonium spectrum under conservation of the spectral 
weight integrated over all charmonia states. In the second class are quark rearrangement (string- 
flip) processes, as indicated in the right panel of Fig. [5] They induce transitions to open charm 
hadrons responsible for charmonium dissociation in hadronic matter, cf. Sect. 21 
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Figure 4: Left panel: The dynamically screened interaction potential V^i^{uj,q) (wavy line), 
determined by the bare potential (dashed line) and the polarization function Ilab{(^,q)- Right 
panel: Cluster expansion for the generalized RPA, when besides free particles (RPA) also two- 
particle states (cluster- RPA) contribute to the polarizability of the medium, see Rcf. [19]. 




Figure 5: Contributions to the dynamical self-energy of a two-particle system in a correlated 
medium at 2"'' Born order. Left panel: impact by two-particle states without constituent 
exchange (van-der-Waals or dipole-dipole interaction). Right panel: constituent-rearrangement 
collisions (string- flip process), from Ref. [53], see also [55]. 



Example 2: String-flip model of charmonium dissociation. Here we give a second 
example for the use of insights from plasma physics by discussing charmonium dissociation 
within the string- flip model of quark matter [T] [S] [271 [2H] ■ In this model the string- type color 
interactions between quarks get saturated within the sphere of nearest neighbors so that in a 
dense system of overlapping quark-antiquark pairs frequent string-flip processes take place in 
order to assure the system is at any time in its minimal energy configuration, see the left panel 
of Fig. E 

When considering a heavy quark-antiquark pair in dense matter with negligible heavy-flavor 
fraction, the Pauli blocking and exchange self-energy contributions are negligible, but the strong 
correlations with light quarks of complementary color within the nearest neighbor sphere will 
result in a meanfield sclfenergy shift (Hartree shift A^) for all quarks [2^. This determines 
the shift of the continuum edge, see the graph (b) in the right panel of Fig. [SI Because of the 
compensation in the Bethe-Salpeter kernel between the effects of the screening of the interaction 
and the self energy shifts calculated with it (see discussion above) , it is suggested that to lowest 
order the bound state energies remain unshifted when increasing the temperature and/or density 
of the medium. In contrast to the first example of Debye screening of long-range Coulombic 
interactions, the screening mechanism in the string flip model is color saturation within nearest 
neighbors, applicable for strong, short-range interactions as appropriate for the case of the 
sQGP at RHIC or dense systems at FAIR CBM. The resulting two-particle energy spectrum 
for charmonium and bottomonium states is shown in the right panel of Fig. [6] where the 
static screening picture (graph (a) as a function of the screening parameter /i = /^^(r) in the 
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Figure 6: Left panel: String-flip process in a dense quark-antiquark system, see also [8| I26|. 
Right panel: Two-particle energies of charmonia and bottomonia states in a statically screened 
potential (a) and in the string-flip model (b), from Ref. [3T] . 

screened Cornell potential [3D] is compared to the string-flip picture (graph (b) as a function 
of the temperature T) , from Ref. [521 ISI] • From this Figure one can read-off the in-medium 
lowering of the dissociation threshold fcg which is the energy difference between the considered 
bound state level and the continuum edge shown as the border of the hatched region. This 
lowering of fep""^ with increasing density and/or temperature leads to a strong increase in the 
quarkonium breakup cross sections by thermal impact [31] and to the bound state dissociation, 
even before the binding energies vanish at the critical Mott densities and temperatures for the 
corresponding states. 

4 Charmonium dissociation in a resonance gas 

An interesting phenomcnological guideline for the present discussion is provided by Ref. [32j 
where the authors show that universal J/i/i-hadron breakup cross sections with correct kinematic 
thresholds but otherwise constant at 3 mb for meson and 5 nib for baryon impact are sufficient 
to explain the observed anomalous J/t/i suppression pattern of the NA50 experiment within 
a multi-component hadron gas. For a recent update, see [33| . The question arises: How do 
these assumptions relate to microscopic calculations of charmonium dissociation reactions in 
hadronic matter? Those can generally be divided into two categories, based on hadronic or 
quark degrees of freedom. We will primarily review and compare the studies of processes in a 
mcsonic medium (predominantly composed of pions and rho mesons) within both approaches, 
including a discussion of in-medium effects. 

Historically, a first calculation of the quark-exchange reaction, J/ij] + ti D + D* + c.c, 
was performed in a nonrelativistic quark model |34| based on applications of the diagrammatic 
technique developed by Barnes and Swanson |35| for meson- meson scattering, see also Ref. |36| . 
This calculation showed a strong energy dependence of the cross section with a peak value 
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Figure 7: Born diagrams for quark-exchange processes (denoted as CI, C2 for "capture" and Tl, 
T2 for "transfer" diagrams) contributing to heavy quarkonium dissociation by meson impact 
(left panels) and the corresponding cross section for J /if) breakup by pion impact (right panel); 
from Refs. |371[3S]. Note that in the total cross section also the processes with charge conjugated 
final states are accounted for. 

of about 6 mb at threshold and a fast decrease due to a momentum mismatch in the overlap 
integrals between the meson wave functions. The thermally averaged cross section, which is 
the relevant quantity for estimating the J/tp dissociation rate, was below 1 mb, roughly in 
accordance with phcnomcnological expectations based on the observed suppression in heavy- 
ion reactions at the SPS. These calculations were extended to other light mesons and excited 
charmonia in Refs. |371 138j. where also more realistic quark-interaction potentials have been 
used. In Fig. [7]wc show the diagrams for quark-exchange processes at first Born order, which 
arc classified as "capture" (C) and "transfer" (T) diagrams depending on whether the quark 
interaction can be absorbed into the external meson lines; the latter are to be understood as 
a resummation of ladder-type quark-antiquark interactions. There are cancellations among 
the contributions of the different diagrams due to the small color dipolc of the charmonium 
state. These cancellations reduce the peak value of the cross section to about 1 mb, as also 
illustrated in Fig. [T] An open question in this approach is whether the double nature of the 
pion - a Goldstone boson of the broken chiral symmetry and a strongly bound quark-antiquark 
- has a strong influence on these results. In the nonrclativistic approach the pion emerges due 
to a large hyperfinc splitting in the Fermi-Breit Hamiltonian (as opposed to, e.g., instanton- 
induced interactions), which is not a robust interpretation. Another question concerns the 
applicability of the truncation of the transfer diagram contributions at the first Born order. 
Ladder-type resummations would lead to s- and t-channel I?-meson exchange processes, which 
arc disregarded in the nonrclativistic quark exchange. Finally, these models are beset with the 
"post-prior" problem due to the ambiguity in the ordering of quark exchange and interaction 
lines. 

These problems can be resolved within relativistic quark models developed on the ba- 
sis of Dyson-Schwinger equations |39j for applications to the charmonium dissociation prob- 
lem j40l 1411 142j . In this approach, the meson-meson interactions arc represented by quark- loop 
diagrams with three (triangle) and four (box) meson legs. The appearance of meson-exchange 
processes can be understood as a ladder resummation of quark interaction diagrams in s- and 
i-channels, see Fig. [8] The results for the J/ if) dissociation cross section by pion impact within 
the relativistic quark model [41], shown in Fig.[8l basically confirm those of the nonrclativistic 
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Figure 8: Left panel: box and triangle diagrams for meson-meson interaction vertices contribut- 
ing to J breakup by meson impact in the relativistic quark model [40| (upper part), and the 
origin of £)-meson exchange from ladder-type resummation (lower part). Right panel: cross 
section for pion induced J /if) dissociation into open-charm mesons (dotted line) composed of 
subprocesses with different final-state D-meson pairs: D + D (solid), D + D* (dashed), D* + D* 
(dash-dotted); from Ref. [5T]. The total cross section includes also the subprocess with D + D* 
final state for which the cross section is identical to the charge conjugated one (dashed). 



approaches [Ml Ell [38] shown in Fig. [3 residual differences may be traced back to the treat- 
ment of the transfer- type diagrams. In the relativistic treatment these diagrams are resummed 
beyond the first Born approximation and result in £)-meson exchange diagrams which are ab- 
sent in the nonrelativistic models. This difference may explain small differences in cross section 
results from both approaches. 

Already in 1998, Matinyan and Miiller [43j pioneered an alternative approach to charmonium 
absorption by light mesons on the basis of an effective meson Lagrangian, with a local ?7(4) 
flavor symmetry strongly broken by the pseudoscalar and vector meson mass matrices. This 
initial version of the chiral Lagrangian approach gave rather small cross sections, cr^^, « 0.3 mb 
at threshold. It turned out [JH [JS] that triple vector- meson couplings and contact terms, 
not included in Ref. [43], result in an increase of the breakup cross section by up two orders of 
magnitude, with a rising energy dependence. A problem for the chiral Lagrangian approaches is 
the treatment of hadrons as pointlike objects, which, as usual for an effective theory, becomes 
unreliable at high momentum transfers (due to quark-exchange substructure effects). The 
composite nature of hadrons can be accounted for in an approximate way by vertex form 
factors [lHl[lS]j which reduce the magnitude, and affect the energy dependence, of the cross 
sections, see Fig. [9] The choice of the cutoff parameters in the form factors has a large effect 
on the charmonium breakup cross sections and remains a matter of debate [48l [49l [50l [51] . 
Progress may be made by calibrating the form factors with microscopic approaches as, e.g., the 
nonrelativistic or the relativistic quark models |52j . While the calculations with the relativistic 
quark model are very cumbersome, the chiral Lagrangian models offer a very effective tool 
to assess many other dissociation processes required for phenomenology, once the formfactor 
question is settled. In particular, these applications include dissociation processes by nucleon 
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Figure 9: Feynman diagrams for J /ip breakup by meson impact in the chiral Lagrangian ap- 
proach (left) and the resulting cross sections for pion- and rho- meson induced processes (right), 
from Ref . [52] . The left panel illustrates the interrelation of the different approaches to charmo- 
nium breakup cross sections: The nonrelativistic potential model diagrams (a) can be redrawn 
as quark loop diagrams with Born order insertions (b). The latter can be either absorbed into 
the nonperturbative meson-quark-antiquark vertices (c, upper line) or after partial resumma- 
tion of all ladder-type diagrams of the Born series redrawn as t-channcl and s-channel meson 
exchanges (c, lower lines). The locan limit of the rclativistic quark model diagrams (c) leads to 
the chiral Lagrangian diagrams (d) with formfactors originating from the quark loop diagrams 
for the meson vertices in (c). 



impact [53|, I54j and bottomonium dissociation |55j . 

The energy dependent cross sections for heavy quarkonia dissociation by hadron impact 
enable to evaluate the temperature- (and density-) dependent dissociation rates in hadronic 
matter. Assuming the (short-distance) vertex functions not to be altered by the surround- 
ing medium, there remains the issue of mass and widths changes of open-charm hadrons with 
temperature and density. These, in particular, imply modifications of the thresholds for the 
breakup processes [Ml [59l [60l [611 |62l [H^ . In Fig. [10] we show the diagrammatic representation 
of the quark-exchange contribution to the J /ip self-energy which develops an imaginary part 
(determining its width or inverse lifetime) due to the coupling to open-charm mesons. For 
the interaction vertex C/ox three different approaches have been discussed above and the corre- 
sponding vacuum cross sections are shown in Figs. [7][9] The theoretical basis for the discussion 
of quark exchange effects to the self-energy of heavy quarkonia in strongly correlated quark 
matter comes from systematic cluster expansion techniques developed in the context of plasma 
physics. For details, see the next subsection and Figs. [2] [3] and [5] Fig. [TT] illustrates the effect 
of the spectral broadening of D-mcsons at the chiral phase transition due to the opening of 
the decay channel into their quark constituents (Mott effect) for temperatures exceeding the 
D-meson Mott temperature T^°^^ k, 172 MeV. Due to an effective lowering of the J /ip breakup 
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Figure 10: Left panel: Diagrammatic representation of the quark exchange process contribution 
to the heavy quarkonium self-energy in a mesonic medium [M]; right panel: same as left panel 
but for quark exchange with baryons in the medium |57] . Three different approaches to the 
interaction vertex Ucx are discussed in the text and shown in Figs.[7][9l Diagrams of this type 
appear in the cluster expansion for two-particle properties, see Fig. O 




. 14 *^ 0.1 0.15 

I 2 T [MeV] 



Figure 11: Left panel: energy (s) and temperature (T) dependence of the effective cross sec- 
tion (a*) for J/ip breakup by p-meson impact. We display p'^a*{s;T) for better visibility of 
the effective lowering of the breakup threshold when temperatures exceed the D-meson Mott 
temperature T^°** w 172 MeV; right panel: temperature dependence of the thermally averaged 
J/4' breakup cross section in a ir-p meson gas; the calculation with vacuum Z?-mesons (dashed 
line) is compared to one with an in-medium broadening of the D-meson spectral function (due 
to the Mott effect at the chiral phase transition) which exhibits a steplike enhancement (solid 
line) caused by the effective lowering of the breakup threshold; from Ref. [56] . 



threshold the temperature dependence of the thermally averaged J ftp breakup cross section in a 
TT-p meson gas exhibits a steplike enhancement |56j . This effect has been discussed as a possible 
mechanism underlying the threshold-like anomalous suppression pattern of J/ip^s observed by 
the CERN NA50 experiment [61] [64j and should also play a role in explaining the final state 
interactions of heavy quarkonia produced in the RHIC experiments. 
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5 Perspectives 



In this lecture we have adapted a general thermodynamic Green function formalism developed 
in the context of nonrelativistic plasma physics for the case of heavy quarkonia states in strongly 
correlated quark matter. Besides the traditional folklore explanation of eharmonium suppresion 
by (Debye-) screening of the strong interaction, we discuss further effects of relevance when 
heavy quarkonia states propagate in a medium where strong correlations persist in the form of 
hadronic resonances. These effects may be absorbed in the definition of a plasma Hamiltonian, 
which was the main result of this contribution. This plasma Hamiltonian governs the in-medium 
modification of the bound state energy levels as well as the lowering of the continuum edge which 
leads not only to the traditional Mott effect for the dissociation of bound states in a plasma, but 
can also be applied to calculate the in-medium modification of quarkonium dissociation rates in 
a consistent way. A detailed recent review |65] summarizes the phenomenological applications 
to heavy quarkonia production heavy- ion collision experiments at CERN and RHIC. Further 
developments of the presented approach shall include in particular applications to quarkonia 
production in dense baryonic matter such as envisaged for the PANDA and CBM experiments 
at FAIR Darmstadt and possibly also for NICA at JINR Dubna. 
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